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Abstract 

Radiatively-driven transfer flow perpendicular to a luminous disk is examined in the subrelativistic 
regime of (v/c) 1 , taking into account the gravity of the central object. The flow is assumed to be vertical, 
and the gas pressure is ignored, while internal heating is assumed to be proportional to the gas density. 
The basic equations were numerically solved as a function of the optical depth, and the flow velocity, the 
height, the radiative flux, and the radiation pressure were obtained for a given radius, an initial optical 
depth, and initial conditions at the flow base (disk "inside"), whereas the mass- loss rate was determined as 
an eigenvalue of the boundary condition at the flow top (disk "surface"). For sufficiently luminous cases, 
the flow resembles the case without gravity. For less-luminous cases, however, the flow velocity decreases, 
and the flow would be impossible due to the existence of gravity in the case that the radiative flux is 
sufficiently small. Application to a supercritical accretion disk with mass loss is briefly discussed. 
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1. Introduction 

Mass outflow from a luminous disk is a clue to the for- 
mation mechanism of astrophysical jets and winds in the 
active objects. In particular, in a supercritical accretion 
disk, the mass-accretion rate highly exceeds the critical 
rate, the disk local luminosity exceeds the Eddington one, 
and mass loss from a disk surface driven by radiation pres- 
sure would take place (see Kato et al. 1998 for a review 
of accretion disks). 

So far, radiatively driven outflow from a luminous 
disk has been extensively studied by many researchers 
(Bisnovatyi-Kogan, Blinnikov 1977; Katz 1980; Icke 1980; 
Melia, Konigl 1989; Misra, Melia 1993; Tajima, Fukue 
1996, 1998; Watarai, Fukue 1999; Hirai, Fukue 2001; 
Fukue ct al. 2001; Orihara, Fukue 2003), as on-axis jets 
(Icke 1989; Sikora et al. 1996; Renaud, Henri 1998; Luo, 
Protheroe 1999; Fukue 2005a), as outflows confined by a 
gaseous torus (Lynden-Bell 1978; Davidson, McCray 1980; 
Sikora, Wilson 1981; Fukue 1982), or as jets confined by 
the outer flow or corona (Sol et al. 1989; Marcowith et al. 
1995; Fukue 1999), and by numerical simulations (Eggum 
et al. 1985, 1988). Line-driven outflows were also exten- 
sively examined (Drew, Verbunt 1985; Vitello, Shlosman 
1988; Proga et al. 1998; Murray, Chiang 1995; Proga et 
al. 2000; Proga, Kallman 2004). In almost all of these 
studies, however, the luminous disk was treated as an ex- 
ternal radiation source, and the radiation transfer in the 
flow was not solved. 

Radiation transfer in the disk, on the other hand, was 
investigated in relation to the structure of a static disk 
atmosphere and the spectral energy density from the disk 
surface (e.g., Meyer, Meycr-Hofmcister 1982; Cannizzo, 
Wheeler 1984; Shaviv, Wehrse 1986; Adam et al. 1988; 



Hubeny 1990; Ross ct al. 1992; Artemova et al. 1996; 
Hubeny, Hubeny 1997, 1998; Hubeny ct al. 2000, 2001; 
Davis et al. 2005; Hui et al. 2005; see also Mineshige, 
Wood 1990). In these studies, however, the vertical move- 
ment and the mass loss were not considered. 

Recently, mass outflow as well as radiation transfer have 
been examined for the first time in the subrelativistic and 
fully relativistic cases (Fukue 2005b, c). In these stud- 
ies, the gravity of the central object and the gas pressure 
were both ignored as a first step. In this paper, we thus 
consider the radiatively driven vertical outflow - moving 
photosphere - in a luminous flat disk within the frame- 
work of radiation transfer in the subrelativistic regime of 
(v/c) 1 , while taking into account the gravity of the central 
object, although the gas pressure is ignored. 

In the next section we describe the basic equations in 
the vertical direction. In section 3 we show our numerical 
examination of the radiative flow. In section 4 we briefly 
apply the present model to the case of a supercritial ac- 
cretion disk. The final section is devoted to concluding- 
remarks. 

2. Basic Equations 

Let us suppose a luminous flat disk, inside of which 
gravitational or nuclear energy is released via viscous heat- 
ing or other processes. The radiation energy is trans- 
ported in the vertical direction, and the disk gas, itself, 
also moves in the vertical direction due to the action of ra- 
diation pressure. For simplicity, in the present paper, the 
radiation field is considered to be sufficiently intense that 
the gas pressure can be ignored (cf. Fukue 2005b): tenu- 
ous cold normal plasmas in the supcr-Eddington disk, cold 
pair plasmas in the sub-Eddington disk, or dusty plasmas 
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in the sub-Eddington disk. As for the order of the flow 
velocity v, we consider the subrelativistic regime, where 
terms of the first order of (v/c) arc retained, in order to 
take into account radiation drag. Under these assump- 
tions, the radiation hydrodynamic equations for steady 
vertical (z) flows without rotation are described as fol- 
lows (Hsieh, Spiegel 1976; Kato et al. 1998 for the full set 
of basic equations) . 

The continuity equation is 

pv = J (= const), (1) 

where p is the gas density, v the vertical velocity, and J 
the mass-loss rate per unit area. The equation of motion 
is 

dv GMz 



K abs 



dz 



(R- 



>R 



[F-{E + P)v], (2) 



where M is the mass of the central object, R = \/r 2 + z 2 , 
r being the radius; r g (= 2GM/c 2 ) is the Schwarzschild 
radius, K a bs and K sca are the absorption and scattering 
opacities (gray), E is the radiation energy density, F is the 
radiative flux, and P is the radiation pressure, measured 
in the inertial frame, respectively. As for the gravity, we 
adopt the pseudo-Newtonian potential (Paczyhski, Wiita 
1980). 

When the gas pressure is ignored, the advection terms 
of the energy equation are dropped (cf. Kato et al. 1998), 
and heating is balanced with the cooling, 



= 9 -p J - cn ahs E + K abs - 



2Fv 



(3) 



where q + is the internal heating and j is the emissivity. 
The last term in the parentheses is the work done by the 
radiation field under the present approximation. 

Under the assumption of steady one-dimensional flow, 
and under the present approximation, the radiation field 
equations (Kato et al. 1998) become 



dF 

dz 



j CK absE + (fvabs " 



Fv 



(4) 



dP pv 



CK a bsE)-p 



^abs 



dz C* " C 

and the closure relation becomes 



[F-(E + P)v},(5) 



P= l -E- 
3 



AFv 
3~^' 
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Eliminating j and E with the help of equations (3) and 
(6), and introducing the optical depth by 



dr = —{k 



abs 



i)pdz, 



(7) 



we can rearrange the basic equations up to the order of 

{v/c) 1 as 



T dv 
CJ Tr 

dz 

' ~d7 

dF 
dr 



GMi 



K abs 



1 



(R-r g ) 2 R 



— (F — 4Pv) , 



K abs 



(tabs " 



F V -, 
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(8) 
(9) 
(10) 



dP 



F - APv 



( K abs 4" K Bca )p C 

In this paper we assume that the heating q + 
portional to the gas density p for simplicity (cf. 
2005b), and set the following relation: 

q + Fo 

1 

TO 



[Vol. , 

(11) 

is pro- 
Fukuc 

(12) 



i^abs ~r 

where Fq is a parameter that determines the magnitude 
of heating, and therefore the magnitude of the radiation 
field. 

We solved equations (8)-(ll) for appropriate boundary 
conditions. 

As for the boundary conditions, we imposed the follow- 
ing cases. At the flow base (disk "inside") with radius r 
and an arbitrary optical depth To, the flow velocity is zero 
(v = 0), the height is zero (z = 0), the radiative flux is zero 
(F = 0), and the radiation pressure is Po, where subscript 
denotes the values at the flow base. At the flow top (disk 
"surface"), where the optical depth is zero, the radiation 
fields should approximately have the values above the lu- 
minous infinite disk, and therefore satisfy the condition 
cP s /F s = 2/3, where subscript s denotes the values at the 
flow top (Fukue 2005b). Hence, the parameters are r, To, 
Fq, and Pq at the flow base, and the mass-loss rate J is 
determined as an eigenvalue of the boundary condition at 
the flow top. 

In numerical calculations, we start from the asymptotic 
solutions near to tq: 



t = T - T, 
1 ^0 

v ~ 

2cJ tq 



F< 



6(fC a bs + K sca )c J 2 T 
^0 



V, 



-t. 
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(13) 
(14) 

(15) 

(16) 
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and solve equations (8)-(ll) down to the optical depth to 
t = 0. We then obtain the final values of v s , z s , F s , and P s 
at the flow top. Generally, the final values do not satisfy 
the boundary condition of cP s /F s = 2/3. We thus iterate 
the value of J so as to satisfy the boundary condition. 
Technically, in some cases, we iterate Po instead of J. 

Finally, we normalize the physical quantities in terms 
of the speed of light c, the Schwarzschild radius r g , and 
the Eddington luminosity Le [= ^ircGM / (K a h s + K S ca)L 
since wc consider the luminous disk around a black 
hole. That is, the normalized quantities are v = v/c, 
z = z/r g , F = F/L E /(4Trr 2 g ), P = cP/Le/{*ki%), and 
J = c 2 J / L-e / (ATrr 2 ) . Moreover, the basic equations (8)- 

(18) 
(19) 



(11) are expressed as 



dv 
J Tt = 

dz 
J d^ = 



(R- 
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dF 
~d7 
dP 



Fo 

Tq 



■Fv, 



(20) 
(21) 



F 

= F-4Pv- —v, 
dr t 

where the symbol "hat" (say, v) is dropped. 
3. Radiative Flow under Gravity 



In this section we show the radiative vertical flow in the 
luminous disk under the influence of the gravity of the 
central object. In order to obtain the solution, as already 
stated, we numerically solve equations (18)-(21), starting 
from asymptotic solutions near to r = tq with appropriate 
initial conditions, down to r = so as to satisfy appropri- 
ate boundary conditions there. The parameters are the 
initial radius r on the disk, the initial optical depth To, 
which relates to the disk surface density, the heating rate 
coefficient Fq, which roughly expresses the radiative flux 
at the flow top and is the measure of the strength of ra- 
diation field to gravity, and the initial radiation pressure 
Pq at the disk base, which connects with the radiation 
pressure gradient in the vertical direction and relates to 
the disk internal structure. The mass-loss rate J is deter- 
mined as an eigenvalue of the boundary condition at the 
flow top. 

Several examples of numerical calculations are shown in 
figures 1 and 2. 

In figure 1 we show the final velocity v s (solid curves), 
the flow height z s (solid ones), the radiative flux F s 
(dashed ones), the radiation pressure P s (dotted ones) at 
the flow top, and the mass- loss rate J (chain-dotted ones), 
as a function of Po at the flow base for r = 3, To = 10, and 
several values of Fq. The quantities are normalized in 
units of c, r g , and L~e / (Airr^) . 

As can be seen in figure 1, the mass- loss rate increases 
as the initial radiation pressure increases, while the flow 
final speed increases as the initial radiation pressure and 
the loaded mass decrease. These properties are similar to 
those in a previous study without gravity (Fukue 2005b) . 
In particular, when the disk luminosity is high (figure lb), 
the results arc quite similar to the previous case without 
gravity. In the case where the disk luminosity is low (figure 
la), on the other hand, the final velocity becomes small, 
as easily expected. Furthermore, if the disk luminosity is 
sufficiently small, the radiative flow becomes impossible 
due to the effect of gravity. 

In addition, as can be seen in figure 1, in order for the 
flow to exist, the radiation pressure Po at the flow base 
is restricted in some range similar to the previous case. 
That is, in the previous study without gravity (Fukue 
2005b), the initial pressure Pq is proved to be restricted in 
the range of 2/3 < cPq/F s < 2/3 + To /2 in the case where 
heating is included. In the present case, when the disk 
luminosity (characterized by Fq) is high (figure lb), Pq is 
restricted in the similar range. When the disk luminosity 
is low (figure la), on the other hand, the values of Pq be- 
come small, but is yet restricted in some range. At the 
upper limit of these ranges, the loaded mass diverges and 
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Fig. 1. Final velocity v B (solid curves), flow height z B (solid 
ones), radiative flux F s (dashed ones), radiation pressure 
P s (dotted ones) at the flow top, and the mass-loss rate J 
(chain-dotted ones), as a function of Pq at the flow base. 
The parameters are r = 3, tq = 10, and (a) Fq = 0.1 and (b) 
Fo = 1. The quantities are normalized in units of c, r g , and 
FE/(47rr|). The leftside scale of the ordinates is for quantities 
except for J, while the rightside is for J. 

both the flow final speed and flow height become zero. 
On the other hand, at the lower limit, the loaded mass 
becomes zero and the flow final speed becomes the maxi- 
mum. 

In figure 2 we show the flow velocity v (solid curves), 
the flow height z (chain-dotted ones), the radiative flux F 
(dashed ones), and the radiation pressure P (dotted ones), 
as a function of the optical depth r for r — 3, To = 10, and 
several values of Pq in a few cases of Fo- The quantities 
are normalized in units of c, r g , and Lb/ '(47rr|). 

When the initial radiation pressure Pq at the flow base is 
large, the pressure gradient between the flow base and the 
flow top is also large. As a result, the loaded mass J also 
becomes large, but the flow speed v is small. When the 
initial radiation pressure Pq is small, on the other hand, 
the pressure gradient becomes small, and the loaded mass 
also is small, but the flow speed becomes large. 

When the disk luminosity (characterized by Fq) is large 
(figure 2b), the flow velocity as well as the radiation fields 
become large. When the disk luminosity is small (figure 
2a), on the other hand, the flow velocity and the radiation 
fields are small. 

In the limit of small Fq (see figure 2a), the radiation 
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Fig. 2. Flow velocity v (solid curves), flow height z 
(chain-dotted ones), radiative flux F (dashed ones), and ra- 
diation pressure P (dotted ones), as a function of the optical 
depth t. The parameters are r = 3, to = 10, and (a) Fq = 0.1 
and (b) Fq = 1. In addition, the values of Pq are (a) 0.2840, 
0.3705, 0.4346, and 0.5480 and (b) 1, 2, 3, 4, and 5 from top 
to bottom of v and z and from bottom to top of F and P. The 
quantities are normalized in units of c, r g , and Le / (47rr? ) . 



fields approach the non-relativistic case of v <C c. In such 
a case, equations (20) and (21) are analytically integrated 
to give 



— (TO' 
TO 



cP = cPo 



-)• 

IF 

2Vo {T - 



Tof 



(22) 
(23) 



where F s = Fq in this limit. In this limit, the relation 
derived in the previous study (Fukue 2005b) again holds: 

cPo 2 1 

F s 

Hence, 



= 3 + 2 T °- 



(24) 



rp4 



cir 



1 

2 T0 



(25) 



where T c is the temperature at the flow base (i.e., Po = 
aT^/3), and T e s the effective temperature at the flow top 



(i.e., F B -- 
imation. 



crT cff ). This is just the modified Milne approx- 



4. Radiative Flow in the Supercritical Disk 

In this section we apply the present model to the mass 
outflow in the luminous supercritical accretion disks. 

When the mass-accretion rate M in the disk around 
a central object of mass M highly exceeds the critical 
rate M cr ;t, defined by M C rit = ^e/c 2 , the disk is believed 
to be in the supercritial regime, and the disk luminos- 
ity exceeds the Eddington one. Such a supercritical ac- 
cretion disk, a so-called slim disk, has been extensively 
studied, both numerically and analytically (Abramowicz 
et al. 1988; Eggum et al. 1988; Szuszkiewicz ct al. 1996; 
Beloborodov 1998; Watarai, Fukue 1999; Watarai et al. 
2000; Mineshige et al. 2000; Fukue 2000; Kitabatake et al. 
2002; Ohsuga et al. 2002, 2003; Watarai, Mineshige 2003; 
Fukue 2004) . It was found that the optically-thick super- 
critical disk is roughly expressed by a self-similar model 
(e.g., Watarai, Fukue 1999; Fukue 2000; Kitabatake et al. 
2002; Fukue 2004). Except for the case of Fukue (2004), 
however, many of these analytical models as well as nu- 
merical ones did not consider the mass outflow from the 
disk surface. Hence, in this paper we adopt the model de- 
velopted by Fukue (2004), as a background supercritical 
disk model. 

In the critical accretion disk constructed by Fukue 
(2004), the mass-accretion rate is assumed to be regulated 
just at the critical rate with the help of wind mass-loss. 
Outside some critical radius, the disk is in a radiation- 
pressure dominated standard state, while inside the criti- 
cal radius the disk is in a critical state, where the excess 
mass is expelled by wind and the accretion rate is kept to 
be just at the critical rate at any radius. Here, the critical 
radius is derived as 

9V3a T , r 



167rm p c 



-M 



input 



1.95mr„ 



where Mi nput is the accretion rate at the outer edge of the 



disk, and m = M 



input 



Outside 



the accretion 



rate is constant, while, inside r c 
vary as 

IQ-Kcm-r. 



the accretion rate would 



M(r) = 



Mj npu t ■ 



(27) 



9v / 3o't 

In such a critical accretion disk, the disk thickness H is 
conical as 



H _ 



m 

1+ 20 



(28) 



where C3 is some numerical coefficient determined by the 
similar procedure in Narayan and Yi (1994) for optically- 
thin advection-dominated disks. The second equality 
of this relation comes from the numerical calculation 
(Watarai et al. 2000). Although the mass loss was not 
considered in Watarai et al. (2000), we adopted this re- 
lation as some measure: when the normalized accretion 
rate m is 100, the coefficient ^/c 7 becomes 0.448. 

Furthermore, in Fukue (2004), several alternatives are 
discussed, and some of them gives the physical quantities 
of the critical accretion disk with mass loss as 
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Fig. 3. Several quantities for each radius r: The disk height 
H (dashed curve) and the optical depth r (chain-dotted one) 
are from the critical model, while the height z B and velocity 
n s (solid curves) at the flow top are the results of the present 
numerical calculations. The quantities are normalized in units 
of c and r g . The parameters of the critical disk is fa = 100 
and a = 1. 
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where a is the viscous parameter. In addition, the mass- 
loss rate J per unit surface area becomes 



/ - 2pH 1 dM - 
2nr dr 



Afinput 



1 



(32) 



In the present non-dimensional unit in terms of c, r g , 
and the Eddington luminosity Le, these physical quanti- 
ties are expressed as 



TO = 


16^/6 r- 

Vr, 

a 


(33) 


F s = 


3 \ 


(34) 


Pa = 




(35) 


J = 


1 


(36) 



where the symbol "hat" (say, f ) is dropped. 

Using these relations, we can solve the basic equations, 
and obtain numerical solutions at each radius r. The ex- 
ample in the case of m = 100 and a = 1 is shown in figures 
3 and 4. 

In figure 3 we show several quantities for each radius r: 
The disk height H (dashed curve) and the optical depth 
To (chain-dotted one) are from the critical accretion disk 
model, while the height z s and velocity v s (solid curves) 
at the flow top are the results of numerical calculations. 
The quantities are normalized in units of c and r g . 

As can be seen in figure 3, the flow height z s , where the 
optical depth is zero, is about 1/4 of the disk scale height 
H. Moreover, the final speed v s at the flow top is quite 
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Fig. 4. Flow velocity v (solid curves) and optical depth r 
(dashed ones), as a function of height z for several values of 
r. The values of r are, from left to right, 2, 3, 4, 5, 6, 7, and 
8. The quantities are normalized in units of c and r g . The 
parameters of the critical disk is fa = 100 and a = 1. 

small in this example; v s ~ 0.001. It is noted, however, this 
speed is not the final wind speed from the disk, but the ini- 
tial starting speed of the mass loss from the disk surface. 
With this starting speed, mass loss from the disk surface 
takes place, and is further accelerated in an optically thin 
regime by the intense disk radiation fields. It should be 
noted that the final speed v s gradually increases with de- 
creasing r since the flux F s increases with decreasing r. 
In the innermost region, however, it decreases due to the 
relativistic gravity of the pseudo-Newtonian potential. 

In figure 4 we show the flow velocity v (solid curves) 
and optical depth r (dashed ones), as a function of height 
z for several values of r. The quantities are normalized in 
units of c and r g . 

As can be seen in figure 4, the flow velocity v roughly 
varies as 

vocz 1/2 . (37) 

This behavior is easily understood, since equation (2) be- 
comes 

dv K abs 4 



—F (38) 
dz c 

in the non-relativistic limit without gravity, and easily 
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integrated to give 



1 

— v 

2 



-Fz, 



(39) 



when the flux is constant. 
5. Concluding Remarks 

In this paper we have examined the radiative trans- 
fer flow in a luminous disk in the subrelativistic regime 
of (u/c) 1 , while taking account of gravity of the central 
object. The flow is assumed to be vertical, and the gas 
pressure is ignored for simplicity, while internal heating 
is assumed to be proportional to the density. The basic 
equations are numerically solved as a function of the op- 
tical depth r, and the flow velocity v, the height z, the 
radiative flux F, and the radiation pressure P are obtained 
for a given radius r, the initial optical depth To, and the 
initial coditions at the flow base (disk "inside"), whereas 
the mass-loss rate J is determined as an eigenvalue of the 
boundary condition at the flow top (disk "surface"). For 
sufficiently luminous cases, the flow resembles the case 
without gravity. For less-luminous cases, however, the 
flow velocity decreases, and flow would be impossible due 
to the existence of gravity. 

The radiative flow investigaed in the present paper must 
be a quite fundamental problem for accretion-disk physics 
and astrophysical jet formation, although the present pa- 
per is only the second step (cf. Fukue 2005b), and there 
are many simplifications at the present stage. 

For example, we have ignored the gas pressure. In gen- 
eral cases, where the gas pressure is considered, there usu- 
ally appears sonic points (e.g., Fukue 2002), and the flow is 
accelerated from subsonic to supersonic. In this paper we 
consider a purely vertical flow, and the cross section of the 
flow is constant. If the cross section of the flow increases 
along the flow, the flow properties such as a transonic 
nature would be influenced. 

Fully relativistic effects are also important. In such a 
fully relativistic case, where the flow speed is of the order 
of the speed of light, the boundary condition at the flow 
top should be carefully treated (cf. Fukue 2005c). The 
radiation drag becomes much more important. 

There remain many problems to be solved. 

This work has been supported in part by a Grant-in- Aid 
for the Scientific Research (15540235 J.F.) of the Ministry 
of Education, Culture, Sports, Science and Technology. 
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